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Abstract 

Throughout the history of the study of turbulence in fluid dynamics, there has yet 
to arise a unique definition or theoretical criterion for this important phenomenon. 
There have been interesting conjectures made by Ruelle [2], Muriel [3], and Getreuer, 
Albano and Muriel [6], however, attempting to provide the sufficient criteria for the 
onset of turbulence. In this paper, a classic equation in fluid dynamics, Burger's 
equation, is solved in one and two dimensions, and these conjectures are illustrated. 
This illustration supports these conjectures by showing that the proposed criteria 
do arise mathematically from the solutions of an equation modelling fluid flows. 
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To this day, there is no unique, universally acceptable theoretical criterion for 
the occurrence of turbulence pp. However, Ruelle [2] conjectured that when 
the time derivative of the velocity inside a fluid becomes infinite, the system 
becomes turbulent. A less stringent criterion was proposed by Muriel [3J, which 
is that if the velocity reverses (even slowly), turbulence could result. The latter 
criterion seems to have been observed in pipe experiments [4l5] . Recently, an 
intriguing suggestion has been made by Getreuer, Albano and Muriel [6] that 
a multiple-valued velocity field may represent turbulence. In this Letter, we 
explore each of these conjectures using the example of Burger's equation. 

Burger's equation is a fundamental PDE from fluid mechanics. It can be ob- 
tained from the more general Navier-Stokes equation under the proper set of 
assumptions. It has been extensively studied over the years, both generally 
(as in [7]), and in particular with respect to turbulence [8)19] . 
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In this paper, the stationary states of Burger's equation in one and two di- 
mensions are studied. Upon application of periodic boundary conditions, the 
solutions are shown to become multi-valued, which allows for the possibility of 
velocity reversal. We also obtain the time- dependent solutions, and show infi- 
nite time derivatives are possible, with the proper choice of arbitrary solution 
constants. 

The development of the above results follows. 

The Navier-Stokes equation can be written as follows: 



Here p represents the mass density of the fluid, v is the velocity, p is the 
pressure, T is the stress tensor, and f represents the sum of any other forces. 

We make the following simplifying assumptions: 



In other words, we consider a Newtonian, incompressible fluid ([2]) under no 
external forces ([3]) or pressure gradient (J3J. 

Setting the velocity constant in time 



is equivalent to solving for the stationary states of the system. In Cartesian 
coordinates and considering the 1-D case, we get 





V • T = V 2 v 
f = 0, 
Vp = 0. 



(2) 
(3) 
(4) 




(6) 



a solution of which is given by 
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Imposing periodic boundary conditions 

v(x) = v(x + L) 
yields the value of the constant C\. 



Cx - (9) 



where k is any nonzero integer. Upon substitution, solution (J7J) becomes 



2{ikir ( knx\ 



Here, C 2 is dropped, equivalent to setting v(0) = 0. 

Note that we now have a family of solutions indexed by the integer k. This 
means that the velocity field can now be considered a multi- valued function. 

Considering the same problem in two dimensions, we have the following set of 
equations for the stationary states: 



/ du x du x \ (d 2 u x d 2 u x \ 
P \ Ux dx +Uy dy J ^{dx 2 + dy 2 )' 1 ] 



where the total velocity function is given by u = u x i + u y j. Solutions to this 
system of differential equations are given by 
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u x {x, y)=Ci + C 5 tanh (d + C 2 x + C 3 y) , (13) 

/ \ C 2 C 4 pC 2 C 5 + 2fi (Cf + Cf) , 

y) = — - - tanh(a 1 + C 2 x + C 3 y). (14) 

Setting the arbitrary constants Ci, C2 and C3 to be imaginary, (i.e., Cj = zi^-, 
where i^- is real) leads to the hyperbolic tangent becoming a tangent function: 



u x (x, y) = C A + C 5 tan (K x + K 2 x + iT 3 y) , (15) 
K 2 C A P K 2 C 5 -2p( 
K 3 pK 3 



%(a;, y) = tanh(i^i + K 2 x + K 3 y). (16) 



Boundary conditions along both the x and y directions can now be imposed, 
similar to what was done in the ID case, to obtain 



u x (x,y) = K 4 + K 5 tan^K 1 + ^x+ } ^y ) j, (17) 

LyipH^ +2/i7r((f ) 2 + (-^) 2 )) kn kn 

u y (x, y) = y y L x »\SlJ \LjJl tan(Kl + ^x+ "-fy) 

pKy L X Ly 

_ k x L y K A 
L x hy 



These 2D solutions are now indexed by two integers, k x and k y , and can again 
be considered multi-valued. 

These ID and 2D solutions may then be interpreted to be possible states of the 
system, with the possibility of transitions between them, following Getreuer, 
Albano and Muriel [6j. This interpretation then allows for an explanation of 
velocity reversal (a signature of turbulence) - the system simply jumps from 
a positive solution to a negative one. 

Finally, in comparison to our results, there is yet another multi-valued char- 
acter in one of the branches of the GAM solution. In one branch, which cor- 
responds to our integer labels, a bifurcation of the solution is observed. Such 
a bifurcation is a characteristic signature of chaos, but which we do not find 
in our example. In this example, the multi-valued nature of the solutions is 
not an intrinsic characteristic, but rather due to the imposition of periodic 
boundary conditions. 
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The time-dependent 1-D Burger's equation takes the form 



f dv dv\ d 2 v . . 

p \Yt +v ^ (19) 



with the corresponding solution 



v( x , t) = - ^ tanh(d + C 2 x + C 3 t). (20) 

L-2 P 



Note that with the choice of parameters C\ = 0, C 2 = ikn/L, and C3 = 0, this 
time-dependent solution ( 1201) reduces to the previously obtained result (ITU]) . 
In addition, by making the arbitrary constants Ci, C 2 , and C3 imaginary, as 
before, we obtain an expression involving the tangent instead of the hyperbolic 
tangent: 



v(x,t) = -^-^^ta,n(K 1 + K 2 x + K 3 t). (21) 
K 2 p 



The time derivative of this solution is 



dv 2pK 2 K 3 



3 (, , ^„„2 



^ p 



1 + tan^x + K 2 x + tf 3 t) , (22 



which evidently approaches infinity for particular values of x and t. 
In 2D, the time-dependent Burger's equations are 

(du x du x du x \ (d 2 u x d 2 u x \ 

P \^ +U ^ +U ^)=»\-M + W)' (23) 

9 \dt +Ux dx +Uy dy ) * \dx 2 + dy 2 ) ' 1 ] 



with the corresponding solutions 
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u x (x, y) = _ PCzC G + 2^Ct + Cj) (c , + ^ + + 

u v (x, y) = C 5 + C 6 tanh (d + C 2 x + C 3 |/ + C A t) . (26) 



We may again choose imaginary constants within the argument of the hyper- 
bolic tangent to transform the function into the tangent. The time derivatives 
of the transformed solutions are then 



— — = 1 1 + tan (A i + K 2 x + K 3 y + K 4 t) 

at pK 2 v ; 

(27) 

= K 4 K 6 (l + tan 2 {K x + K 2 x + K 3 y + K A tj) . (28) 



These time derivatives again clearly become infinite at certain values of x, y, 
and t. 

The points at which the time derivatives of these 1-D and 2-D solutions become 
infinite may very well be points where the onset of turbulence occurs, as 
hypothesized by Ruelle [2]. 

In conclusion, we exhibit an example of the Ruelle conjecture, illustrate the 
modification proposed by Muriel, and quite interestingly, satisfy the conjecture 
of Getreuer, Albano and Muriel that multi- valued velocity fields could possibly 
describe turbulence. 
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